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ABSTRACT 


This paper discusses a procedure by which the absolute alignment of the 
magnetic axes of a triaxiai magnetometer sensor with respect to an external, 
fixed, reference coordinate system, can be accurately determined. 

The procedure does not require that the magnetic field vector 
orientation, as generated by a triaxiai calibration coil system, be known to 
better than a few degrees from its true position, and minimizes the number of 
positions through which a sensor assembly must be rotated to obtain a solution. 

Computer simulations have shown that accuracies of better than 0.4 
seconds of arc can be achieved under typical test conditions associated with 
existing magnetic test facilities. 

Tne basic approach Is similar In nature to that presented by MePherron 
and Snare 097B) except that only three sensor positions are required and the 
system of equations to be solved is considerably simplified. Applications of 
the method to the case of the MAGSAT Vector Magnetometer are presented and ti e 
problems encountered discussed. 


INTRODUCTION 

The problem of determining the absolute orientation of a magnetic field 
vector has been solved traditionally by assuming that the field orientation 
can be accurately established by the geometry of a calibration coil. This 
method Is generally sufficient to determine sensor orientations to within a 
few minutes of arc from Its true direction but if higher accuracies are 
required not only roust we take into account additional parameters In the coil 
gecroetry and Its construction, but also \ts time and temperature stability. 

A straightforward method of determining the absolute orientation of a 
magnetic field vector is by rotation of a fluxgate magnetometer sensor (or any 
other vector sensor) on a surface about an axis approximately parallel to the 
field direction. The surface orientation can then be adjusted to obtain a 
constant reading. Independent of rotation angle. It can easily be shown tliat 


litider these conditions the surface nonnal, as measured with respect to a 
reference coojlnate system, is parallel to the applied magnetic field vector. 
It ir obvious that for maximun seo'^itivity, the field must be applied 
approximately normal to the sensor axis. 

A major drawback of tMs approach Is the fact that accurate planar 
rotations can only be obtained about the local vertical axis due to gravity 
effects on the sup^wt structures. In addition the method is extremely lime 
consuming and requires complex sup^v^rt instr urentat ion such as non-magnet ic , 
preclsioti, j-degrees-of-f reedixn fixtures of pe:'formance Ci^mparable to first 
order theodo L i tes . 

The pi'oocdurc presented beli'w onviates the need for multiple rotations 
and special mounting futures and allows the stmul tatieous determinat ion of the 
sensor assembly and test coil facility al igtiment para u*ters by means of a 
simple Iterative algorltUn and mevisurt'ments obtained for three discrete and 
fixed sensor positions. Variations with time andd tero|.x'rature of test coll 
orientations do not affect Uic detorm inat ion of the sensor alignment if their 
character Istic time is long comjvired with the test time. Thus significant 
savings in test time are realised since many s»^urces of rror which otherwise 
would have to be accounted for in a rittu'r com pi 'x a!id iruer acting scheme, arc 
el iminat ed , 

ai.k;mmknv wr^.KMl’^A^d^v MKrn.'.p 

The basic assumpt iot^ in the metlu'^i is t iiat A' the deviations fr^'m 
or thogona I i I y ot' t.iu' scniaor assembly .ind to*t v'o i I s .v'c :^al I (within a few 
degrees> the :neasurt>metU s ot^tained frvxn a :v'nsor ;uount<\l nonnal to the field 
direction will refleot- th<‘ sum ot t ne dev i ■; ions of uv' sensor and test coil 
ay st tMn . 

A "‘et'eretice coi’^rdinate syst e:n is first aconrately established by a pair 
of fif'st order t iu'Ovio i 1 1 es whioii are rigidly mounted to suitable suppi^rts , 
levi^lcvi and *'ofer'enoed to a si aL' I d'ocotioji, Ttie coil sy:^ m 

orientation is approx imat e i y al ignel witti the referc'iicc ooordlnate systecn and 
ttie magnetic field veotors ;n W\o reference and c'vl coor'dinate systems are 


laled by 


( 1 ) 


where [B] = [I] (unit matrix) if the systems are nearly aligned. Thus, in 
general, « 1 for i j. 

In analogous fashion we define a matrix CA] which relates the 
measurements in the sensor coordinate system to the reference coordinate system 

Mg = [A] Mj^ (2) 

♦ 

where Mg is the measured field vector. Again, if we assune near alignment, 
lAijI << 1 for i ^ j and 1. 

The problem then reduces to determining the elements of [A] from a set of 
measurements obtained by varying il and [A] in a known way. This can be 
accomplished by energizing one axis of the coil system at a time and 
reorienting the sensor assembly to exchange rows or columns of [A]. 


in general then we have 


fig = [A] = [A] [B] 


(3) 


Since the coils are energized one axis at a time and restricting ourselves to 
unit magnitude vectors, we can write as a unit matrix 




c 



1 0 0 
0 1 0 
0 0 1 


[I] 


( 4 ) 


It then follows that the measurements can be organized in a 3 * 3 matrix such 
that 


[Mg] = [A] lB] (5) 

Under the assumption that the off-diagonal elements of [A] and [B] are small, 
we can rewrite (5) as 
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V W / 


[Mg] =■ LI + 6A] [I ♦ dBj 
and neglecting second order terras. 


CMoJ » Cl] [6AJ > [SB] 


U) 


It l3 worthwhile to express ij) in corop<'>nent form and we shall use the 
uppersv-ript (1) to indicate that these values correspond to the first sensor 
position (axes nearly parallel to test coils axes) 
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The 


first subscript of denotes the sensor axis being read while the second 
denotes the coil system axis ( i .e . , x, y, z) which is energized. The same 
convention applies to the elements of the matrices [A] and [B], 


If we now rotate the sensor assembly about the z-axis of the reference 
system, by exactly 90**, the matrix [A] will now take the form 


lA j 


. (J) 


"^xy 

"*xx 


'^y 

‘^yx 

^yz 

~^zy 

’.x 

^ZZ 

u 



tiie resultant 

matr ix . 


(d' 


(?) 

wnere we nave denoted uy lAj 
oan be easily aooomplished by ineans of reference theodolites and optical cubes 
mounted on the sensor assembly, as described later on in this paper. 


If we now energize the coil system axes in tne same order as before we 


have 





( 10 ) 


[Bj 


or, in component form 


[M.J 
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(-a ♦b ) (a ♦b ) 
xy yx xy yx 


‘V^y’ 


(^a ♦b ) 
xy XX 
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‘"zz'^zy* 


(n) 


Where the upperscript (^) denotes sensor position number 2. Equations 19) and 
(11) constitute a system of 12 equations in 12 unknowns but as shown by 
McPherron and Snare (1978) the characteristic matrix of the system is 
singular. Thus an additional sensor rotation is required. We choose to 
rotate the sensor exactly 90“ about the X-axis of the reference system. Hie 
sensor alignment matrix then becomes 


[A] 


O) . 


-8- 
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xy 

XX 

XX 
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yy 
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zy 

XX 

XX 


( 12 ) 


Energizing the coil system axes in sequence we obtain 


tMgJ 


( 3 ) 


(3) 

(3) 

<i)i 1 
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^x 

”xy 

«xz 




(3) 

(3) 
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^y 



-1 

^”*yz“^xy^ ^®yx^^xz^ 

(3) 

(3) 
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”zy 
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‘-‘zy-^yz’ 



(13) 


where the upperscript (j) denotes sensor fv>sitlon number 3. It is clear that 
equations (8), (ID and (13) allow the twelve coefficients a^j, b^^j (i W j) to 
be estimated. The diagonal elements a^^ and b^^ have to satisfy the direction 
cosine constraints. 
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(14) 


a,. = [1 




“u = " 




and hence can be calculated from the solution of (d) , ( 11 ), and ( 13 ) 


Note that it is not necessary to simultaneously solve these equations, 

since id) and ( 11 ) allow a , a , b , a , a , b and b to be determined 

iy £y A xz y^ 

While the remaining coefficients can he determined using (11) and (13). 


It is then convenient to express the solutions as a set of four systems 
of linear equations, as follows 

1-1 0 1 1 1 ■' 


CCi^J 


- 1 0 0 1 


0-1 0-1 


0-110 


10 10 


~ ( 1 ) 

M 




10 0 1 


0-1 1 0 


0 1 


\ 1 
zy ( 


and 



^23 ^24 ®33 


(C^4 


t1 - 

[1 - (C^^ ^ C^^)] 


These matrices [A] ^ and [B] represent our first attempt to determine the 
exact solution and hence, at this stage of the calculation we will have 

(i) 

i IA]„^ (2 a: 

i = 1,2,3 


To improve our estimates we can utilize (24) to implement an iterative 
scheme which will converge to the desired result ^®^est' 


To accomplish this we first form the estimated matrices 


—AAA 
12 11 13 


-A,, A 


21 23 


' -A A A 

^^32 33 


-^2 

-^3 

^1 

-a,3 

"^23 

^21 

-A . 

-^3 



We then obtain the products 


where the A. .‘s are elements of the matrix [A] . , 

Aj est 


(1) 

^^S^est = 

(1) 

tA^est 

^®^est 


(2) 

^^S^est " 

(2) 

^®^est 

(27) 

^^S^est " 

(3) 

^‘B3est 



and the differences with the measured values 


> 

II 

- 

(1) 

^^S^est 


4^2 = 


(2) 

^^S^est 

(28) 

= 

- 

(3) 
o est 



We then solve equations (14) through (17) utilizing the updated values of 
and defined by 


[Mg] 

(1) 

. /I) 

^^S^est 

♦ 

Afi^ 


[Mg] 

(2) . 

(2) 

^^S^est 

♦ 

AMg 

(29) 

[Mg] 

(3) 

(3) 

^^S^est 


AM^ 


kth iteration 

'‘^-’>th 

iteration 



until the smallest element of , ^ 1^2 or does not exceed a predetermined 
small value typically chosen as 10”^. When this condition is satisfied we 
consider and to be the desired solution. 

Note that in equation (28) the value of [fi and 

o o S 

remain constant durir^ the iteration procedure since they represent the actual 
measurements. 


It is instructive at this p^^int to compare the present method with that 
given by McPnerron and Snare (197 <3). The basic difference is the method of 
solution of the approximate equation (?)• In their paper they state that a 
tuinimun of four different orientativons of tne sensor are required which lead 
to a systesn of 36 equations iti 12 unknowns. The solution of this system is 
iinpiernentec usiOig tne singuisr value deco.npi'^sit ion metnod of Lanzcos. In the 
present method, the solution of equations (15) through (19) is trivial and 
only three seru'or positions are required , considerably simplifying tne 
measuranents task. 

COMPUrEK 


The procedure presented in the preceeding section was programmed on an 
IBM 5100 computer utilizing the APL language (.IBM* 1977). This simulation 
program is snown in Appendix A. Tne inputs to the program were two matrices 
tAJ and [Bj representing the coil and sensor misal ignraents , Representatives 
values for [A] and [BJ were chosen and the theoretical measurement value 
computed. These values were in turn used to recover the alignment matrices 
[A] and [Bj with the metnod presented. 

Two examples are shown in Tables I and II. The first represents typical 
values expected in the case of the MAG3AT Vector Magnetometer and a typical 
coil system. [AJ and [Bj are the input matrices and CSV] and [BV] tne 
corresponding estimated matrices. is can be observed tne input matrices were 
recovered exact to the Ulth decimaj. place in only j iterations. The second 
example 1 , fable II) represents an extrene case where angular misal ignments as 
large as 16.7^ were allowed for tne sensor matrix. Again, the procedure 
recovered the input matrices witii an accuracy 01 ^ x id in on ly 0 

iterations . 


For these two examples we have assumed that ine measurements are error 
free. Beic^w we iiscuss several sources ''f error atid their effect upon the 
overall accuracv of tne alignment determination , 



tHtiUH suuncfca AKU IMtiH iW*ACl UH UtftKALL AtUUKACil 


Two fundamental measurement limitations must be taken into account to 
estimate the overall accuracy of any aligrment method and they are: a) 
resolution and b; noise. We will discuss later additional sources of error 
which fortunately can be adequately accounted for by the alignment 
determination method. 

The resolution with which we can obtain the measurements establishes a 
minimun value of signal that can be reliably detected. For the MAGSAT case 
this was not a limitation because the measurements were obtained with 6 digit 
resolution. Thus in the absence of noise, we theoretically resolved 0.1 uV of 
signal which corresponds to 1 part in 4000 (the calibration constant for the 
instruoent was 4 mv/nT or 2.5 * 10 nT in a field of 50,000 nT, which 
corresponds to an angular error of 1 x 10~^ seconds of arc. However, noise , 
in both the coil system and magnetometer, constitutes a more fundamental 
limitation. Typical values of noise were 0.1 nT RMS for the coil system 
"zero" fi^ld and 0.01 nT RMS for the magnetometer. Thus, it became extremely 
important to minimize the calibration facility noise. 

For purposes of discussion let's assune that the "noise" (including all 
contributions) can be reduced to 0.05 nT by suitable procedures with a worst 
case value of 0.1 nT. If tne test field is 50,000 nT, this implies that 
angular deviations smaller than 0.2 seconds of arc, with a worst case value of 
0.4 seconds of arc, cannot be reliably detected. 

These effects were simulated in the computer program where artificial 
random "noise" in multiples of 0.4 arc second peak amplitude w^s added to the 
theoretical measurements. The results of these simulations are presented in 
Tables III and IV for the same cases previously presented in Tables I and II. 
The nunber of iterations in this case was fixed at six since obviously the 
10”^ bound for could not be achieved. As expected, the alignment matrices 
obtained show deviations of the same order of magnitude as the random noise 
amplitude. 
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An additional source of error is tne accuracy with which the test field 
can be established. Fortunately a proton precession magnetometer used to 
determine the absolute value of the field to better than ] nT which, for the 
small misalignment angles involved, does not affect the values obtained to any 
significant extent. Ht^wever, to illustrate its contr ibution , let us assume 
that the sensor and coil system are misaligned by O.i® and the field is in 
error by ^ nT, For small angles 

la a (30) 

where a r misalignment angle 

Aa = deviation of a from true value 


z applied field 


AH 


c 


deviation of applied field from true value. 


The s^rrective angular error inroduceu is then 0.04 arc seconds for a 
b0,000 nT test field. As indicated by (30) this error is propc^r t ional to the 
misalignment angle ana hence they should be made initially as small as 
possible. Inis also implies tnat the axis tne reference coordinate system 
should be aligned with tnr coil system aXiS as accurately as possible. 


An additional source of error wtiich must be minimized is tnat due to 
field gradients associated witC; tne coii system. Since tne sensor assembly 
generally cannot rotate about the sensor axes, each rotation results in a 
si ight translation for at least ^ o: i};rec sensors constituting tne 

triad. Tills translation ,..s unimportant unless tne field gradients are 
relatively large in whi ‘h case the sensors are exposed to different fields 
depending on tneir position wittiin the coll system. it is clear from the 
preceeding discussion that the field gradient should be less than 0,1 nT 
a tae largest- dimensicr. sen -or assembly to maintain an overall 

accuracy comparshle 'o tnat Ic: : :: y the noise cha»-'ec ter ist ics of the 


CO ; I system , 



B 


Finally, we must consider the effects introduced by the presence of 
magnetic materials in the Immediate vicinity of the coil system which, by 
induction effects, alter the direction of the *’free space" field produced by 
the coil system. Examples of this problem were the thedolites themselves 
which incorporate in their construction soft magnetic materials with high 
effective permeability. Normally, these instrunents are mounted remote from 
the center of the coll system but unfortunately close enough to introduce 
deviations in the orientation of the generated field of the order of a few arc 
seconds . 

The magnitude of the Induction field can easily be determined with n 
proton precesion magnetometer. Since the alignment determination method does 
not require precise knowledge of the coil system alignment, the theodolites 
can be considered as integral parts of the system since their geometry and 
orientation remain fixed during ail the magnetic tests. This obviously can be 
accomplished with great accuracy since first order instrunents are generally 
used for calibration. 

a pplication to the magsat vector magnetometer 

The sensor assembly for the KAGSAT Vector Magnetometer is shown in Figure 
1* Three ring core fluxgate sensors are mounted orthogonally on a glass- 
ceramic base. Two optical cubes were bonded to the sensor assembly as 
illustrated in the figure and the relative alignment of the cube faces 
measured with a set of first order theodolites. A description of the 
magnetometer overall design and electronics has been given by Acuna et al . , 
1976 and Acuna, I960 so it will not be repeated here. 

The calibration and alignment tests were conducted at the Goddard Space 
Fiight Center 6 meter Magnetic Test Facility. Two concrete piers separated by 
approximately 400 meters established an optical azimuth reference, as shown in 
Figure 2. This base'ine was used to verify the absolute azimuthal orientation 
of two first order theodolites mounted inside the coil building which 
established the primary, orthogonal, reference coordinate system. The system 
was verified periodically by means of auxiliary mirrors and found to yield 
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repeatable results with a typical uncertainty of less than 2 arc-seconds. The 
facility was equipped with a state-of-the-art servo control system to monitor 
and correct earth's field variations in the range of 0 to 25 Hz. The typical 
noise level of the system was 0.1 nT RMS. All calibration activities took 
place during selected "quiet" evenings to minimize errors due to refraction 
effects by hot air currents and other significant disturbances such as nearby 
electric railroad traffic and wind induced distortions of the coil building. 
The computer programs described in Appendix 'A' were modified and expanded as 
shown in Appendix 'B' for these tests. 

During engineering model tests the effects of cross-field non- 1 .nearities 
in ring-core fluxgate sensors were evaluated and found to be significant for 
fields greater than 1b, 000 nT. This problem is shown in P'igure j. Ring core 
fluxgate magnetometers with Large effective ( t/d) ratios exhibit deviations 
from linearity of the order of 1 part in 10 fields applied parallel to 

the sensitive axis. This figure degrades significantly wtien a large field is 
simult.- leously applitnl in a direction transverse to the sensing axis, as shown 
in Figure j. This is due to the appearance of large amplitude second harmonic 
signals at the sensor terminals which are in phastr quadrature with the signal 
produced by the oti-axis field. These large signals affect the linear 
operation of the electronics and lead to the observed instriment response . 

Nv'te that the respt^nse function depends u|.x'n the applied on-axis and 
cross-ax is fields. 

The strategy followed for MAdSAT was to mathemat ical i y model out these 
non-1 inear ities rather than correct them in the itistrimentat ion . To determine 
the dependence of the alignmerjt angles with tne amplitude of the test field 
caused by the non-linear re 3 p«^nse, the method descriceo in the previous 
section was used with tests fields of ^y.),000 nV. 3‘>.000 nT and Vj,000 nT. As 
expected, the results obtained for tne sensor aligiment matrix did vary with 
the amplitude of the test field, but not those for the test coil system 
aligiment matrix. Ttus, of course, is wtiat would be anticipated since the 
coil system alignment matrix is test- field amplitude independent. 

One impAMaant fact derived fi'ixn the measurements was that only those 
matrix elements associated with dirocixons in tne plane of a given sensor ring 





core were field anplltude dependent. This is not surprising since in the 
direction perpendicul jr to the plane of the core there exists a large 
deaagnetizing factor due to the narrow ring-core- sensor geometry and resulting 
small (l/d) ratio (Bozorth, 1951; Acuna, 1969). 


Thus we can write 


[A] = 


XX 




xz 


*yx^®x*^y^ *yy 


yz 


zx 




zz 


(31) 


The functions A^j constitute second order corrections to the basic 

measurements and hence are not strong functions of B. Thus for all practical 
purpo-»es choosing Bj^ = ^ (raw) in ( 31 ) above does not introduce any 

significant errors in the determination of the A^j terms. Typical results 
obtained for A^z's as a function of the test field amplitude are given in 
Table 5 . 


Tne measurements were least squares fitted to functions of the form 


A . , a . » B . 

A (B. ,B.) = Sin 

^ B. A. 

J ij 


(32) 


for A^ 2 * ^21 ^31 excellent results, liow (32) only models the 

response to cross-axis fields without regard for the magnitude of the on axis 
field. Hence (32) must be expanded to include this dependence 



The functional form of was determined experimentally by full M t 

steradian mapping of the instrunent response function. The following 
functional relation was found to fit the measurements with the required 
accuracy 


», (B.) . C * 0 » — ^ . C, exp(C B — J-) (34) 

‘ 0 ' 1 iBJ ' i ‘ IBjl 

The final values derived for the HAGSAT flight sensor coeffi- ants are given 
in Table 6. The final absi’luLe alignment accuracy achieved through the 
measurement and modeling activity was estimated at t 3 arc seconds. 

SUMMARY 


A relatively simple procedure to accurately determine the absolute 
alignment of vector magnetometer sensors has been presented. The method 
minimizes the nunber of sensor orientations necessary to obtain a solution and 
requires simple mathematical operations. Computer simulations have demon- 
strated the rapid convergence of the solutions to exceedingly small values of 
error, even for initial deviations from orthogonality as large as 16 degrees. 
Several error sources limiting the obtainable accuracy in practical 
applications were presented and it was shown that angular determination 
accuracies of the order of 0.4 arc seconds are technically achievable. 

Finally, the application of the method to the MAGSAT Vector Magnetometer 
alignment was presented including second order effects associated with large 
(t/d) ring core fluxgate sensors. 
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FIGURE CAPTIOHS 


Figure 1 - Schematic representation of the MAGSAT Vector magnetometer triaxlal 
sensor assembly. Three ring-core fluxgate sensorc 're mounted orthogonal 
to each other. Two optical cubes bonded to the assembly define the 
reference coordinate system. 

Figure 2 - The MAGSAT optical reference system as implemented at the NASA-GSFC 
6-meter Magnetic Test Facility. The 400-meter baseline stability was 
checked periodically against stellar references. 

Figure i - Response of large (t/d)-ratio ring-core fluxgate sensors to cross 
and on-axis fields simultaneously, applied in the plane of the core. Tne 
deviations from linearity are produced by large quadrature signals 
generated by the sensor under these conditions (large external fields). 

For external field < 5000 nT the effect is negligible. 


TABLE 1 


MAGSATSIM 0 


ENTER SIMULATED SENSOR ALIGNMENT MATRIX 
U- 

A 

ENTER SIMULATED COIL ALIGNMENT MATRIX 


0 : 

B 

A2Y,BZX> AYZ^BXZ :0. 0089995 ”0 . 000367‘m"0 . 0070028 "0 
AZY.BZX. AYZ.BX2 :0. 0089984 *0 . 00036861 ~0 . 007004 "0 . 
AZY,BZX. AYZ.BXZ ;0. 0089984 ~0 . 00036861 "0 . 007004 *0 . 


A 

0.99999 
0 .0056 
0 .0001 
B 

1 . OOOOEO 
5.3000E~5 
"3.6700E"4 
SV 

9.9999E"! 
5.6000E‘3 
1 . 0000E“4 
BV 

1 .OOOOEO 
5.3000E*5 
•3.6700E‘4 


0 . 002 
0 . 99996 
0.009 

2. 0000E"4 
1. OOOOEO 
3.8900E‘4 

2. OOOOE“3 
9,9996E*^1 
9.0000E'3 

2.0000E'4 
1 . OOOOEO 
3.8900E'4 


"0.003 
"0 . 007 
0.99996 

"1 .2000E"4 
4.5000E"4 
1 . OOOOEO 

"3. OOOOE‘3 
"7. OOOOE‘3 
9.9996E‘l 

"1 .2000E"4 
4.5000E"4 
1 .OOOOEO 



r B 


.00012528 

00012403 

00012403 
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TABLE II 


MAGSATSIM 0 

ENTER SIMULATED SENSOR ALIGNMENT MATRIX 
A 

ENTER SIMULATED COIL ALIGNMENT MATRIX 
B 

AZY,B2X, AYZ.BXZ "0.11 1 . 1 154E "5 "0 . 30 0 03 ~4 
AZY,EZX, AY2-BX2 "0.10997 3 . 9936E "5 "0 . 300 05 
AZY>H2X, AYZ.BXZ : '0.10997 4 . 0 079E "5 "0 . 30 0 05 
AZY.PZX. AYZ.DXZ ."0.10997 4 . 0274E "5 "0 . 30 0 05 
AZY.BZX, AYZ,BXZ ; '0.10997 4 . 0278E ' 5 *0 . 30 0 05 
AZY,B2X, AYZ,BXZ :'0. 10997 4 . 028E "5 "0 . 30 0 05 


A 

0.97468 
0 . 05 
0 . 005 
B 

1 . OOOOEO 
3. 0000E"6 
1 . OOOOE'5 
SV 

0.97468 
0.05 
0 . 005 
BV 

1 . OOOOEO 
2.9993E"6 
1 . 0000E"5 


0 . 1 

0. 95263 

~ 0.11 

5. OOOOE‘4 

1 . f OOEO 
4.5000E~5 

0 . 1 

0.95263 

~ 0 .11 

5. OOOOE~4 
1 . OOOOEO 
4.5000E~5 


"0 . 2 
"0.3 
0.99392 

4.0000E~5 
2. OOOOE~4 
1 .OOOOEO 

" 0.2 
~0 .3 
0.99392 

4 .OOOOE'5 
2. 0000E"4 
1 . OOOOEO 


7675E"6 

~8. 0389E"6 

~8.6579E"6 

~8.7035E~6 

~8.7113E'6 

8.712E~6 
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TABLE III 


MA6SATSIM 1 

ENTER SIMULATED SENSOR ALIGNMENT MATRIX 

U 

A 

ENTER SIMULATED COIL ALIGNMENT MATRIX 


AZY.BZX. AYZ. 
AZY.BZX, AYZ, 
AZY,BZX. AYZ> 
AZY.BZX, AYZ, 
AZY.BZX. AYZ. 
AZY.BZX. AYZ. 
A 

0 . 99999 
0 . 0056 
0.0001 
B 

1 .OOOOEO 
5.3000E'5 
"3.6700E"4 
SV 

0.99999 

0.0055999 

0.0001 

BV 

1 .OOOOEO 
5.2955E*5 
“3.6688E'4 


BXZ : 0.0089995 "0 . 0 0 036744 "0 . 0070028 ' 
BXZ : 0.0089982 ~0 . 00036875~0 . 007004 't 
BXZ -0.008998 "0 . 0C036899"0 . 0070039 "t 
BXZ : 0.0089977 ~0 . 00036924"0 . 0070039 ' 
BXZ ; 0.0089975 '0 . 0 0 036949 "0 . 0 070038 ' 
BXZ : 0.0089972 "0.00036974*0.0070038 ' 


0 .002 

0.99996 

0.009 

2.0000E"4 
1. OOOOEO 
3.8900E‘4 

0.0019997 

0.99996 

0.0090001 

1 .9978E~4 
1 .OOOOEO 
3.8909E*4 


0 . 003 
0 .007 
0.99996 

"1 .2000E‘4 
4.5000E'4 
1 .OOOOEO 

■0.0029998 

"0.007 

0.99996 

"1 .2017E"4 
4.5002E"4 
1 . OOOOEO 


0.00012528 
.00012416 
. 00012412 
0 . 00012407 
J . 00012403 
0 .00012398 


TABLr IV 


MAGSATSIM 1 

ENTER SIMULATED SENSOR ALIGNMENT MATRIX 
U- 

A 

ENTER SIMULATED COIL ALIGNMENT MATRIX 

Q 

B 

A2Y,BZX/ AYZ/BXZ 1 . 1154E~5~0 . 30003 *4.7675E"6 

AZY.BZX, AYZ^BXZ : ~0 . 10997 3 . 999iE 'S-O . 30005 '8.2689E"6 
AZY,BZX, AYZiBXZ •‘0.10997 4 . 0 1 17E ‘5 ‘0 . 30 0 05 ‘9.0327E~6 
AZY,BZX. AYZ.BXZ : ‘0,10997 4 . 0295E “5 ‘0 . 30 0 05 ‘9.2234E‘6 
AZY . B2X, AYZ, BXZ : ‘0 . 10997 4 . 02&4E ‘5 ‘0 . 30 0 05 ‘9.376E‘6 
AZY.BZX, AYZ.BXZ ;‘b.i09?7 4 . 027E "5 ‘0 . 30 0 05 ‘9.5215E‘6 
A 

0.97468 0.1 ‘0.2 

0.05 0.95263 *0.3 

0.005 ‘0.11 0.99392 

B 

l.OOOOEC 5.0000E‘4 4.0000E‘5 

3.0000E 6 l.OOOOEO '20000E~4 

l.ODOOE'5 4.5900E‘5 l.OOOOEO 

SV 

0,97468 O.l *0.2 

0.05 0.95263 ‘0.3 

0.005 ‘0.11 0.99392 

BV 

l.OOOOEO 4.9978E‘4 3.9951E‘5 

3.1344E‘6 l.OOOOEO 2.0003E‘4 

1.0051E‘5 4.5191E‘5 l.OOOOEO 
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Test Field Amplitude [nX] 


Coefficient 

15.000 

35,000 

55,000 

*.y < V» 

- 1.0604 X 10“^ 

- 9.26198 X 10"^ 

- 7.10277 X 10-** 

*xz 

7.92 X 10"^ 

7.92 X 10"** 

7.92 X 10"** 

*,x <V<» 

- 2.868 X 10“^ 

- 2.61318 X 10"^ 

- 2.18877 X 10"^ 

*,x 

2.208 X 10"^ 

2.208 X 10"^ 

2.208 X 10"^ 


2.3325 X 10"^ 

2.3325 X 10"^ 

2.3325 X 10"^ 

A,„ (B,=0) 

- 2.84214 X 10"^ 

- 2.65283 X 10"^ 

- 2.33114 X 10"^ 

zy 2 
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APPENDIX B 




7MEASc:n:i7 

7 A HE AS B,P1,P2,P3 

hihis program assumes ihai both i-ielii pol, are applied eor each se'nsgr posihon 

[2] fi A IS THE TEST FIELD AMPLITUDE AND B ARE THP: CALIBRATION FACTORS T VOLTS/ GAMMA I X.Y.Z 
i:3J C03: 'ENTER MEASURED VALUES POR POSITION ttl ( X -Hi -U ) , 0 , D-U , U-E , N - S , X , Y . Z , ORDER' 

1:4:] PI‘-P1- 7 3 p(Pl*- 7 3 pPl«-n)Il,:i 

1:'=;] 'DIFFERENCE MATRIX FOR POS.«l ';P1 
[.6J 'ARE VALUES OK'^ ENTER YES OR NOP' 

1: 7T 4( (+/□=' NOP' )^^-^3)/C03 

tei NORMALIZED MATRIX FOR POSITION «I ' , ( PI <-NMAG PI) 

i:V] CQ4 ENTER MEASURED VALUES FOR POSITION < X->D U . Y->N - S ) ' 

[10 1 P2<-P2 7 3 p(P2<- 7 3 pP2< [])I. l;T 

i:ilT 'DIFFERENCE MATRIX FOR POSITION «2 ' ; P2 
ri21 ARE VALUES OK'? ENTER YES OR NOP' 
i:i3j ->( (+/n='NOP ‘ )=3)/C04 

114] 'NORMALIZED MATRIX FOR POS , «2 ' , < P2<-i9MAG P2) 

i:i5J C05 'ENTER MEASURED VALUES FOR POS . 113 ( X->N -S - Y-»U-D ) ' 

L16] P3<-P3- 7 3 P<P3<- 7 3 pP3<-f])Ll,] 

i:i7] 'DIFFERENCE MATRIX FOR POS . il3 ' , P3 
[18] ARE VALUES OK"? ENTER YES OR NOP' 

C19] -X (-»/n= 'NOP ■ )=--3)/C05 

C20] NORMALIZED MATRIX FOR POS . «3 ' ; < P3<-NMAG P3) 

1:21] B1<-AP<-P1 

[22] B2«-AP90<-P2 

i:23] B3«-AP90P^P3 

C24] 1 MAGSATSIM 0 

C25] 'SENSOR ALIGNMENT MATRIX' ;SV 
L26] 'COIL ALIGNMENT MATRIX ',BV 
7 


V 

[ 1 :i 
1:23 

C33 

C43 

c:53 

C63 

L73 

[;8:i 

C93 

C103 

i:in 

r.123 

1 : 133 

C143 

ri53 

II163 

C173 

[183 

ri93 

2 0 3 
[: 2 1 3 
1:223 

C233 
1:24 3 
1251 
C263 
C27 3 
T283 
T29 3 
i:;3()3 
[’313 
i:;32 3 


, «• , .3 . «= , H , u , i,. , U3 , K . . ► 

ftA IS A CONTROL. 1 FOR USE WITH HE.AS, 0 FOR STAND ALONE. 

KH- 3 3 P 0 1 1 1 0 1 1 1 0 
K2<- 3 3 .o 0 1 1 1 1 0 1 0 ^ 

K3<- 3 3pll0011101 
«K IS RANDOH NOISE IN MULTIPLES OF 0.4 

ENTER SIMULATED SENSOR ALIGNMENT MAIRIX’ 

A«-D 

ATI ; 13 <-( l-(Ar.l ; 23*2) +AI.1 ;33*2)*0 .S 

AC 2 , 23 «-( 1 -(AL 2 ; 13 * 2 ) +Ai:2, 33 * 2) *0 .U 
A[ 3 . 33 ^( 1 -(AC 3 , 13*2)+A[:3, 2 3 * 2)*0 .S 

■ENTER SIMULATED COIL ALIGNMENT MATRIX 

B<-n 


k 1 ; 1 3<- (1 - ( DC 1 , 23 *2 ) +DT 1 , 3 3*2 ) *0 
BC2 ; 23^(1 - (BC2;1 3*2 )+BC2; 33*2 )*0 
BI:3;33<-<1-(BC3, i:i*2)+B[ 3,23*2)*0 

AP<-A+ . *B 

A90C , 33<--< A90<-®( lOA) ) T , 33 
B2^AP90«-A90+ . *B 

A90PI. ; 1 33<— (A90P< 1<I>A) I , 1 33 

B3<-AP90P<-A90P+. *B 
Bl«- 3 3 p( ,Bi)+KxlE '8x9?2f. 

B2«- 3 3 p( ,B2)+KxlE '8X9?25 
B3< 3 3 p( ,B3)+KxlE'8x9?25 
C02 :M1«-0M1<- 44pl0100l01 
M2<-BM2<- 44pi010010 1100 
M3<-BM3<- 44p0”1100 '10 1 10 


0 1 1 
1 
1 


'10 0 1 
1 1 0 
10 0 1 


M5<-0M5<- 44p010'110 1 001 10 1 00 1 

COl^ VAU APCl . 23 , API.1 ; 33 , AP90C 1 , 211 - AP90T 1 : 33 
VA':’<-APC2 ; 1 3 » APt3 , 1 3 » AP90 [. 2 ; 1 3 » AP90 C 3,13 
VAkAP90t;3,33,AP90P[:3i 1 3 , AP90Pr:3 , 33 . AP90r.3 ; 13 


:2,K3 


c ;i n VA'.i<- Apvof :> , ;i , apvoi i,2i. apvopi; . 21 , Apyopf 1 .. 21 

I: J4J A3.^-M1 •* . XVAI 
C35 3 A2<M2+.xVA2 

i:36) A3<M3+.XVA3 

[373 . XVA5 

C38.1 5V< 3 3 f>l , ,A1L1D. - < (AlC23 + A5r.lD)- 2) , .A2I1 1. < 1, .A5C2D- , ( (A2C23+A3C1D)- 2) , /A3C23, . 1 

[373 BV«- 3 3 pj . , ( (A1C33 + A5C33)^2) , , A1C43. . +( (A2C33 + A3i:‘f3)^-2) . . I . . A3L'33. . A2C+3 . , A5C43. . 1 

I'vO] SVU , 13^"( l-(SVr 1 ,23«2)+SVri ,33«2;*0,5 

[ 413 SVf. 2, 23»- ( 1-<SVC2; I 3»2)+SVC2, 3 3*2)«0 . 5 

r423 SVr;3;33<- <l'(SVr3, n«2)+SV(:3,23*2)»0 .5 

1:433 BVCl, 13< 23»2)+BVr 1,33*2) «0, 5 

C44 3 BV[:2,2J* <1-<BVC2; 1 3*2 ) +BVL2 , 33«2 ) #0 . 5 

L453 BVt3,33< <1--(BVL3, 13*2>+&VC3,23*2)*0 .5 

t463 •A13.B12. A31,B21 ^ASCl 33.A3II1 43 

L4 73 AP«-AP+(Ll<-KlxBl SV+ , xBV) 

C483 D+ ( 3oL1.)x360x360Q- 02 

L49 3 A7()t: , 33*-- (AVO*-«< 14>SV) )C ,33 

C’503 AP90< AP90+(L2<-K2xB2-A90+ . xBV) 

[;S13 n<- ( ‘3 oL2)x3600x360^o2 

C523 A90PI;; 1 33* - < A90 Pf 14»SV ) I , 1 33 

i C;533 AP90P<^ AP90P+(L3<-K3xB3 -A90P+ . XBV) 

1 rr»43 D<-< "3oL3)x3600x360rO2 

f C'353 'SV ,SV,f.lAVCl573, ’BV ,BV 

t T563 UI + 1 

I rS73 ♦(I==7)/0 

I ['383 360 0x360-: 02) X ' 3o I MC 1 19 I , ( ,L1) , ( ,L2) . ,L33>il)/C01 

!■ V 


vMAGrnjv 

V 7<-hAG K 

I n 7^ K 7 i f>h 

r?] 71 < I t'pzr . n > , ni 7i:: ' i ti»zc , 2 j ) , 2]<-zc ( i f^'zc . 33) . 3.i<-a 

Z 33 zt;-; 1 tizi , 11 ), i:i<^zcat 4 >zi: , 2 J) ,23< zc( it^zr , 33 ) , 33 ^ -a 

r '+ '1 2* 3 f> / V 2 p , 6 3 p3i,Z)-2 

t G ] Z' Z 3 3 p .Pi 

V 


Vt: PRO R 1.113 V 

7 A r PROP H , SV? . SV3 . T1 , T2 , 13 

t:i3 HH 15 THt: TEIST F IC LI' AMf'LITUDE AND E< ARt THE CALIBRATION FACTORS 
[ 21 3 3 p< AxBE I T ) , 0 , 0 , 0 . < AxBE 23 ) , 0 , 0 , 0 , AXBC3 1 

L3J ■POSITION «i Rf:C0N5T. MFA5 ' 

T 4 3 n< T 1 < NB t . X 5V+ . x BV ) 

C53 5V2C , 3 3«- - ( S'72< D ( 14)GV ) > [ , 33 

L63 POSITION »2 REICONST. MEAS ’ 

173 ri< T2< NB+ . X (SV2+ , xBV) 

[83 SV3C, 1 33^ (SVT^-KDSV)!, 1 33 
i::?] 'POSITION «3 RECONGT. rtEAS ' 

t;i 0 3 f|< T3<- NB+ . X ( SV3+ . xpV > 
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